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Abstract 



QCD sum rules for decuplet baryon two-point functions are investigated using 
a comprehensive Monte-Carlo based procedure. In this procedure, all uncer- 
tainties in the QCD input parameters are incorporated simultaneously, re- 
sulting in realistic estimates of the uncertainties in the extracted phenomeno- 
logical parameters. Correlations between the QCD input parameters and the 
phenomenological parameters are studied by way of scatter plots. The pre- 
dicted couplings are useful in evaluating matrix elements of decuplet baryons 
in the QCD sum rule approach. They are also used to check a cubic scal- 
ing law between baryon couplings and masses, as recently found by Dey and 
coworkers. The results show a significant reduction in the scaling constant 
and some possible deviations from the cubic law. 
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I. INTRODUCTION 



The QCD sum rule method |L] is a powerful tool in revealing the deep connection between 
hadron phenomenology and QCD vacuum structure via a few condensate parameters - 
vacuum expectation values of QCD local operators. This nonperturbative method has been 
successfully applied to a variety of problems to gain a field-theoretical understanding into 
the structure of hadrons (for a review of the early work, see Ref. 0), and continues to be 
an active field ||. 

The decuplet baryons have been studied in the past using the QCD sum rule method. 
In Ref. H, A was first studied. It was later extended to S* in Ref. 0. In Ref. ||, sum rules 
for the decuplet family were given (they were also reported in Ref. f2|). Ratio method was 
used there to extract the masses. No attempt was made to extract the current couplings. 
No anomalous dimension corrections were considered. In Ref. 0, A was studied including 
leading-order a s corrections. In Ref. ||, A was studied. A perusal of these works reveals 
discrepancies among the QCD sum rules derived where comparisons are possible. The 
analysis methods employed were relatively crude. Often a 10% or better accuracy was 
claimed in all the extracted quantities without the support of rigorous error analysis. 

In this work, we decide to re-derive the sum rules for all members of the decuplet family, 
consistently including operators up to dimension 8, first order strange quark mass correc- 
tions, flavor symmetry breaking of the strange quark condensates, anomalous dimension 
corrections, and factorization violation of the four-quark condensate. Furthermore, we try 
to assess quantitatively the errors in the phenomenological parameters, using a Monte-Carlo 
based procedure 0. This procedure incorporates all uncertainties in the QCD input pa- 
rameters simultaneously, and translates them into uncertainties in the phenomenological 
parameters, with careful regard to OPE convergence and ground state dominance. The goal 
is to get a realistic understanding of the predicative ability of the standard implementation 
of the QCD sum rule approach for the decuplet baryon two-point functions. Of particular 
interest is the baryon coupling to its current. This quantity is crucial to studying matrix 
elements of these baryons, such as magnetic moments, transition moments, axial charges, 
tensor charges, etc., because they utilize the coupling as normalization. 

II. METHOD 

The starting point is the two-point correlation function in the QCD vacuum 

KM =ijd A x e^(0 | T{ % (x) 7j„(0) } | 0), (1) 

where r/ M is the interpolating field (or current) with the quantum numbers of the baryon in 
question. The interpolating field excites the ground state as well as the excited states of 
the baryon from the QCD vacuum. The ability of the interpolating field to annihilate the 
ground state baryon into the QCD vacuum is described by a phenomenological parameter 
\b (called current coupling or pole residue) 
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(0 1 77 J .Bps) = Xbu^p, s), 



(2) 



where w M is the Rarita-Schwinger spin-vector. This parameter plays an important role in 
evaluating matrix elements of the baryons. 

The lowest dimensional interpolating fields for the decuplet are uniquely defined. As- 
suming SU(2) symmetry in the u and d quarks, they can be written as: 



<(x) = e abc (u aT (x)C lfl u b (x)) u c (x), 



(3) 



7% {x) = V /T73e a6c [2 (u aT {x)C^s b {x)) u c (x) + (u aT (x)C^u b {x)) s c (x) 



(4) 



V ? (x) = sjljie abc [2 (s aT {x)C ltl u b {x)) s c (x) + (s aT \x)C lfl s b (x)) u c (x) 



(5) 



rj°(x) = e abc (s aT (x)C lfl s b (xj) s c (x) 



(6) 



Here C is the charge conjugation operator and the superscript T means transpose. 

The QCD sum rules are derived by calculating the correlator in ([!]) using Operator- 
Product-Expansion (OPE), on the one hand, and matching it to a phenomenological repre- 
sentation, on the other. The obtained tensor structure has the form 



n^„(p) = fliO 2 ) + U 2 (p 2 ) g^p + 



(7) 



where the hat notation denotes p = p a j a . Two QCD sum rules can be derived from the two 
invariant functions Ui(p 2 ) and Ui(p 2 ) for each member using standard procedure. The sum 
rule from Ux is called chiral-odd since it involves dimension-odd condensates only. Similarly, 
the sum rule from 112 is called chiral-even since it involves dimension-even condensates only. 
These two structures are considered because they receive contributions from spin-3/2 states 
only. The chiral-odd sum rules at the structure g^ are given for A: 
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E L 2 / 27 M 2 -—ab L 16 / 27 = Ai M A e - M V M \ 



(8) 



for £*: 



\ (f s + 2)aE 1 L 16 / 27 M'- 2 - Us + 2) m 2 a E L 2 / 27 M 2 - i (/. + 2) a b L 16 / 27 
+1 m s E 2 L- 8 / 27 M'+ 2 m s K v a 2 L 16 / 27 = A 2 , e~ M ^' M \ 



(9) 



for S*: 



~ (2f s + l)aE 1 L 16 / 27 M'-\ (2f s + 1) m\a E L 2 / 27 M 2 — ^- (2f s + l)ab L 16 / 27 
9 9 54 

+m s E 2 L" 8 / 27 M 6 + ^ m s f s n v a 2 L 16 / 27 = \l« M.- e - M l*l M \ 



(10) 
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and for f2: 

if.aEk L^ 27 M*- 2 -f s mla E L 2 / 27 M 2 - A f s a b 

+ ^ m s E 2 M 6 + 2 m s / 2 a 2 L 16 / 27 = A 2 M n e"^ 2 . (11) 

The chiral-even sum rules at the structure g^p are given for A: 

- E 2 L 4 / 27 M 6 — —bE L 4 / 27 M 2 + - ^a 2 L 28 / 27 - - mla 2 L u ' 27 ^- = A| e ~ M i/ M2 , (12) 
5 72 3 9 M 2 A ' v ; 

for £*: 

1 £ 2 L 4 / 27 M 6 - A h E L 4 / 27 M 2 + i (2f s + 1) ^a 2 L 28 / 27 - A (2/. + 1) m 2 a 2 L 14 / 27 A, 
+i m s (4 - /.) a £ L 4 / 27 M 2 - A ™ s (14 - 5 /.) m\a L~ 10 / 27 = A 2 , e"^ 2 , (13) 

for E*: 

\ E 2 L 4 / 27 M 6 - A 6Eo L 4/ 27 M2 + 4 /s(/s + 2) ^ L28/27 _ 7_ Ws + 2) m2fl2 Ll4/27 _1_ 
+§ m s (/. + 2) a £ L 4/27 M 2 - 1 m s (2/. + 7) m 2 a L" 10 / 27 = A 2 , e ~ M ^ M \ (14) 
and for f2: 

\ E 2 L 4 / 27 M 6 - A 6 Eo L 4/ 27 M2 + 4 /2 ^ L28/27 _ 7 /2 m2fl2 Ll4/27 _1_ 

+3 m s /. a £ L 4/27 M 2 - | m s /. m 2 a L" 10 / 27 = A 2 e~W . (15) 

In the above equations, a = — (2n) 2 (uu), b = {g 2 G 2 ), (ug c cr ■ Gu) = — m 2 , (ww), \b = 
(2"7t) 2 Ab. The ratio f s = (ss) / (uu) = (sg c cx-Gs) / (ug c a-Gu) accounts for the flavor symmetry 
breaking of the strange quark. The four-quark condensate is not well-known and we use the 
factorization approximation (uuuu) = k v {uu) 2 and investigate its possible violation via the 
parameter k v . The anomalous dimension corrections of the various operators are taken 
into account via the factor L = [a s (/i 2 ) / 'a s (M 2 )} = \\n{M 2 / K 2 QCD ) / \n(pi 2 /A 2 QCD ) I , where 
H = 500 MeV is the renormalization scale and Aqcd is the QCD scale parameter. As usual, 
the excited state contributions are modeled using terms on the OPE side surviving M 2 — > oo 
under the assumption of duality, and are represented by the factors E n (x) = \ — e~ x J2 n x " I m 
with x = w 2 B /M B and % an effective continuum threshold. Note that wb is in principle 
different for different member of the decuplet and we will treat it as a free parameter in the 
analysis. 

These sum rules constitute a complete, up-to-date set of QCD sum rules for the decuplet 
baryon two-point functions under standard implementation of the approach. The reader 
can find differences in a number of Wilson coefficients and in the anomalous dimension 
corrections, when comparing with the existing calculations. 
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III. ANALYSIS 



The basic steps of the Monte-Carlo based analysis are as follows 0. Given the uncer- 
tainties in the QCD input parameters, randomly-selected, Gaussianly-distributed sets are 
generated, from which an uncertainty distribution in the OPE can be constructed. Then a 
X 2 minimization is applied to the sum rule by adjusting the phenomenological fit parameters. 
This is done for each QCD parameter set, resulting in distributions for phenomenological 
fit parameters, from which errors are derived. Usually, 100 such configurations are suffi- 
cient for getting stable results. We generally select 1000 sets which help resolve more subtle 
correlations among the QCD parameters and the phenomenological fit parameters. 

The Borel window over which the two sides of a sum rule are matched is determined by the 
following two criteria: a) OPE convergence — the highest-dimension-operators contribute no 
more than 10% to the QCD side; b) ground-state dominance — excited state contributions 
should not exceed more than 50% of the phenomenological side. The former effectively 
establishes a lower limit, the latter an upper limit. Those sum rules which do not have a 
valid Borel window under these criteria are considered unreliable and therefore discarded. 
The emphasis here is on exploring the QCD parameter space via Monte Carlo. The 10%-50% 
criteria are a reasonable choice that provide a basis for quantitative analysis. Reasonable 
alternatives to the criteria are automatically explored in the Monte-Carlo analysis, as the 
condensate values and the continuum threshold would change in each sample. 

The QCD input parameters and their uncertainty assignments are given as follows. The 
quark condensate in standard notation is taken as a = 0.52 ± 0.05 GeV 3 , corresponding 
to a central value of (uu) = — (236) 3 MeV 3 . For the gluon condensate, early estimates 
from charmonium [l[j place it at b = 0.47 ± 0.2 GeV 4 , a value commonly used in QCD sum 
rule analysis. But more recent investigations support much larger values [p!0| -[T3||. Here we 



adopt b — 1.2 ± 0.6 GeV 4 with 50% uncertainty. The mixed condensate parameter is placed 
at itIq = 0.72 ± 0.08 GeV 2 . For the four-quark condensate, there are claims of significant 
violation of the factorization hypothesis |T0|-p!2]|. Here we use k v — 2 ± 1 and 1 < k v < 4. 
The QCD scale parameter is restricted to Aqc , d=0.15±0.04 GeV. We find variations of 
Aqcd have little effects on the results. The strange quark mass is taken as m s = 0.15 ± 0.02 
GeV. The value of f s has been determined in and is given by f s = 0.83 ± 0.05 after 
converting to our notation by 7 = f s — 1. These uncertainties are assigned conservatively 
and in accord with the state-of-the-art in the literature. While some may argue that some 
values are better known, others may find that the errors are underestimated. In any event, 
one will learn how the uncertainties in the QCD parameters are mapped into uncertainties 
in the phenomenological fit parameters. 

In illustrating how well a sum rule works, we choose to plot the logarithm of the two 
sides of the sum rule against the inverse of M 2 . In this way, the right-hand side will appear 
as a straight line whose slope is — M\ and whose intercept with the y-axis gives a measure 
of the coupling strength. The linearity (or deviation from it) of the left-hand side gives an 
indication of OPE convergence and the quality of the continuum model. The two curves 
should match in the defined Borel region for a good sum rule. Fig. [IJ shows such a plot 
for the chiral-odd sum rules at the structure g M „. The resulting fit parameters are given in 
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TABLE I. Monte-Carlo analysis of the chiral-odd sum rules for the decuplet. The third column 
is the percentage contribution of the excited states to the phenomenological side at the lower end of 
the Borel region (it increases to 50% at the upper end). The results are obtained from consideration 
of 1000 QCD parameter sets. 



Sum Rule 


Region 


Cont 


w 


A 2 


Mass 


Exp. 




(GeV) 


(%) 


(GeV) 


(GeV 6 ) 


(GeV) 


(GeV) 


A 


0.95 to 1.17 


31 


1.65± 0.22 


2.26± 0.89 


1.43± 0.12 


1.232 


S* 


0.82 to 1.29 


8 


1.80 


2.83± 0.32 


1.394± 0.052 


1.384 




0.99 to 1.42 


13 


2.00 


4.32± 0.47 


1.505 ± 0.037 


1.533 





1.05 to 1.59 


8 


2.30 


7.19± 0.75 


1.676 ± 0.031 


1.672 


TABLE II. 


Similar to Table 


|, but for the chiral- 


even sum rules. 


J. lit: ocCOHLl IL)W 111 


cdCll oLLlll 


rule shows the results with the leading 


-order a s corrections included. 






Sum Rule 


Region 


Cont 


IV 


A 2 


Mass 


Exp. 




(GeV) 


(%) 


(GeV) 


(GeV 6 ) 


(GeV) 


(GeV) 


A 


1.15 to 1.56 


10 


2.20 


4.13± 0.65 


1.19± 0.13 


1.232 




1.15 to 1.39 


23 


2.00 


3.79± 0.51 


1.23 ± 0.11 




S* 


1.13 to 1.63 


6 


2.40 


5.63± 0.48 


1.36 ± 0.13 


1.384 




1.13 to 1.47 


15 


2.20 


5.49± 0.41 


1.38 ± 0.11 






1.10 to 1.71 


3 


2.60 


7.71± 0.37 


1.52 ± 0.13 


1.533 




1.10 to 1.56 


8 


2.40 


7.92± 0.43 


1.53 ± 0.10 




n 


1.08 to 1.76 


2 


2.70 


9.16± 0.43 


1.61 ± 0.12 


1.672 




1.08 to 1.61 


5 


2.50 


9.54 ± 0.53 


1.61 ± 0.10 





Table |. The experimental masses are taken from Particle Data Group fl4 |. Fig. || shows a 
similar plot for the chiral-even sum rules. The resulting fit parameters are given in Table O. 

First, let us note that the chiral-odd A sum rule (§) is the only one that allows a three- 
parameter search. In the other sum rules, there is not enough information in the OPE for 
a three-parameter search. What happens numerically is that the search algorithm return a 
continuum threshold that is either zero or smaller than the mass, which is clearly unphysical. 
In order to proceed, we decide to fix the continuum thresholds at certain values and perform 
a two-parameter search on the masses and couplings. By adjusting the continuum thresholds 
to values that reproduce the experimentally known masses, the current couplings X B are left 
as predictions from the self-consistency requirement of the corresponding sum rules. We have 
tried a different approach: fixing the masses at known values and search for the continuum 
thresholds and couplings. This was not successful due to the same reasons as above. It is 
satisfying to observe that the obtained continuum thresholds in Table | roughly coincide with 
the first excited state in each channel from Particle Data Group A(1600), S*(1840), 
2*(1950), fi(2250). 

The predicted A mass lies above the experimental value by less than two standard devi- 
ations. This small overestimation of Ma is a typical result in the QCD sum rule approach. 
For example, in Ref. §] it was found M A ~ 1.37 GeV, A| ~ 2.3 GeV 6 , and w=2.2 GeV; 
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5. 

pi 




0.2 0.4 0.6 0.8 1.0 1.2 

If- 2 (GeV~ 2 ) 



1.4 1.6 



FIG. 1. Monte-Carlo fits of the chiral-odd sum rules at the structure g^ v . Each sum rule is 
searched independently. The solid line corresponds to the ground state contribution, the dotted 
line the rest of the contributions (OPE minus continuum). The error bars are only shown at the 
two ends for clarity. 



C 

Pi 




0.2 0.3 0.4 0.5 0.6 0.7 

M~ z (GeV" 2 ) 



0.8 



0.9 



FIG. 2. Similar to Fig. [j], but for the chiral-even sum rules (without a s corrections) at the 
structure g^p. For better viewing, the curves for each baryon have been shifted downward by 0.2 
units relative to the previous one. 
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and in Ref. § it was found M A ~ 1.36 GeV, A| ~ 1.53 GeV 6 , and io=1.58 GeV. In these 
works, the traditional value of b = 0.47 GeV 4 was used. Using this value in the present 
study yields M A ~ 1.385 GeV, A| ~ 2.00 GeV 6 , and w=lM GeV. The couplings from 
these studies are consistent with each other. It is worth mentioning that both M A and \ 2 B 
as determined from the QCD sum rule method agree with those from other calculations. 



For example, a lattice QCD calculation [15[ gives M A — 1.43 GeV, A^ ~ 2.13 GeV ; and an 
instanton liquid model calculation gives M A — 1.43 GeV, \ 2 B ~ 1.70 GeV 6 . The reason 
for the overestimation of the A mass in QCD sum rule method remains an open question. It 
is likely that the power corrections are large in this sum rule and the OPE is not sufficiently 
convergent. On the phenomenological side, one notes that the continuum contribution in 
this sum rule is the largest among the decuplet family, greater than 30% in the entire Borel 
region used. This also indicates that the continuum model may not be adequate. 

Further examination of Fig. [I] and Fig. reveals that the errors in the chiral-odd sum 
rules are generally smaller than those in the chiral-even sum rules. More importantly, the 
convergence of the OPE in the chiral-odd sum rules are better than that in the chiral-even 
sum rules (as evidenced by the deviation of the dotted line from linearity). Consequently, the 
results from the chiral-odd sum rules are more reliable than those from the chiral-even sum 
rules. This fact is a general feature of baryon two-point functions as recently emphasized in 



Ref. ]17| . It can be traced to the fact that even and odd parity excited states contribute with 



different signs. The contributions of positive- and negative-parity excited states partially 
cancel each other in the chiral-odd sum rules, whereas they add up in the chiral-even sum 
rules. Further evidence for the unreliability of chiral-even sum rules will be discussed below. 

The effects of perturbative corrections to first order in a s were also studied in this work. 
These corrections have been calculated in |7J and are given by 

1 }-\ 1E )^ and (i-^)?' (16) 

to the dimension-three (a) and dimension-five {rn^a) terms, respectively, in the chiral-odd 
sum rules @ to (|TT|); and 

/539 1 \ a s , /113 22 \ a. 



1e -f, and - [t^ + ^e -f, (17) 



V 90 3 ' ) 7T ' V108 3 



7T 



to the dimension-zero (the identity operator) and dimension-six (a 2 ) terms, respectively, in 
the chiral-even sum rules ([12]) to (|l^). Here je — 0.58 is the Euler constant. At the scale 
of about 1 GeV 2 , a s /n ~ 0.12. So the leading-order a s corrections amount to about 5% in 
the chiral-odd sum rules and about 70% in the chiral-even sum rules. Their effects on the 
spectral parameters are given in Table SIS db second entry for each of the chiral-even sum 
rules. We find the effects of a s corrections in the chiral-odd sum rules are small and can 
be safely neglected. They were not listed in Table |. On the other hand, inclusion of a s 
corrections in the chiral-even sum rules leads to a shrinkage of the valid Borel regions and an 
increase of the continuum contributions, signs of deterioration of the sum rules. Since the 
identity operator term is closely tied to the continuum model, it turns out that the effects 
of the a s corrections can be compensated by simply shifting down the continuum threshold 
by about 200 MeV. 
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When considering the uncomfortably large a s corrections in the chiral-even sum rules, 
the possibility of a significant dimension-two power correction from a summation of the 
perturbative series |IB| , |l9fl , coupled with the large uncertainties associated with the four- 
quark condensate, one has to conclude that the QCD side of these sum rules are really poorly 
known. The spectral properties extracted from them are most likely unreliable. Therefore, 
we caution against the use of these chiral-even sum rules in favor of the chiral-odd sum rules. 

Since all the parameters in the Monte-Carlo analysis are correlated, one can study the 
correlations between any two parameters by looking at their scatter plots. Such plots are 
useful in revealing how a particular sum rule resolves the spectral properties. Fig. ^ shows 
the correlations of vacuum condensates with the A coupling for the chiral-odd sum rule (^j). 
Similar plots hold for the mass and the continuum threshold and are not shown. The 
uncertainties in the quark condensate reveal anti-correlations with the spectral parameters, 
while those in the gluon condensate and the mixed condensate display positive correlations. 
The fact that all of the condensates have significant correlations with the phenomenological 
parameters suggests that all three terms in the OPE side of this sum rule play an important 
role in determining the outcome. This may be an indication that the OPE in this sum 
rule is not yet sufficiently convergent. The fact that all three phenomenological parameters 
are correlated with the vacuum condensates in the same manner suggests that attempts 
to fine tune the condensates will increase or decrease the phenomenological parameters 
simultaneously in this sum rule. 





a (GeV 3 ) 



0.0 0.5 1.0 1.5 3.0 2.5 3.0 3.5 



b (GeV 4 ) 




(GeV 2 ) 



FIG. 3. Scatter plots showing correlations between A A and the QCD input parameters for the 
chiral-odd sum rule (||). The results are drawn from 1000 QCD parameter sets. The distributions 
for Ma and w& have qualitatively the same shapes as for A A , and are not shown here. 

Fig. H shows similar scatter plots for Q at the chiral-odd sum rule (|TT|). Here three more 
parameters come into play: k v , m s and f s . It is interesting to observe that the mass and the 
coupling have different correlation patterns. They are opposite for the quark condensate. 
Similar is true for the strange quark mass, although less pronounced. The mass is negatively 
correlated with k v , while the coupling has little correlation with it. The correlations with b, 
ml and f s appear fairly weak. 

We have examined scatter plots for all of the sum rules. We find that qualitatively, 
the strange baryon (£*, S* and Q) sum rules have very similar patterns within the same 
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chirality, despite some subtle differences. The patterns in the chiral-even sum rules are very 
different from those in the chiral-odd sum rules. We will not elaborate further about them. 



O 1.70 
w 1.65 

tfi , „„ 









a (GeV 3 ) 



b (GeV 4 ) m 3 (GeV 2 ) 



m. (GeV) 




FIG. 4. Similar to Fig. 0, but for the chiral-odd f2 sum rule (FT 



In a recent work fl2T| , the authors observed that baryon current couplings and their 
masses obey a simple cubic scaling law: A# ~ Mf , across the octet as well as the decuplet. 
This result is based on the couplings extracted from existing QCDSR calculations. Some 
theoretical justification was given based on general scaling arguments for QCD and a simple 
light-cone constituent quark model. Using our results for the decuplet, we are able to check 
more carefully this claim. 

In Fig. H following Ref. [p0[] , we plot the couplings X% = 2{2n) A \ 2 B as a function of the 
masses, using points extracted by Ref. |H| and our points. The nucleon point is taken from 



Ref. 0] which is the most recent analysis of the nucleon mass sum rule. It lies slightly below 



the point used in Ref. |p0| . For the octet A, £ and 5, there exist no new analyses to our 
knowledge, so the same points are used. The important difference is that our values from 
Table | for the decuplet £*, S* and Q are significantly smaller than the ones extracted by 



Ref. |[20|| . The A points roughly agree. This is true even for the values extracted from the 
less reliable chiral-even sum rules in Table 0, which are slightly larger than those from the 
chiral-odd sum rules. Note that we have taken into account a factor of 2 difference in the 
definition of A^ in the comparisons. We attempted to fit the new points with a cubic law 
and obtained X 2 B = 0.71M| compared to their result of A^ = 1.35M|. Notice, however, that 
there are considerable deviations from the points in the new fit, especially the octet strange 
baryon points. This calls for a re-examination of the QCD sum rules for the octet A, £ 
and H. Given the great importance of a scaling law between baryon current couplings and 
their masses and its phenomenological consequences, more investigations are clearly needed 
to resolve the deviations, perhaps coupled with other methods such as lattice QCD. 
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0.9 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 
Baryon Mass (GeV) 



FIG. 5. Baryon current coupling (X%) versus the baryon mass(M^). The empty squares are 
the points extracted by Ref. pCfl , and the filled circles are from this work. The solid line represents 
the fit to the square points: \ 2 B = (1.16) 2 M^, as obtained in Ref. pQ] . The dashed line is the fit 
to our points: A| = (0.84) 2 M|,. 



IV. CONCLUSION 



We have re-derived and re-analyzed in detail the QCD sum rules for decuplet baryon 
two-point functions using a comprehensive Monte-Carlo based procedure. Predictions were 
obtained for the current couplings with an accuracy on the order of 10% for £*, 5* and Q, 
and 40% for A (see Table |), using realistic error estimates for the QCD input parameters. 
The results are useful in evaluating matrix elements of decuplet baryons in the QCD sum 
rule approach. Our calculations confirmed the general feature that chiral-odd sum rules are 
more reliable than chiral-even sum rules as far as baryon two-point functions are concerned. 
Correlations between the QCD input parameters and the phenomenological parameters are 
studied by way of scatter plots. Some insights are gained on how a particular sum rule 
resolves the spectral properties. Finally, the couplings obtained in the present study are 
used to check a possible cubic scaling law between baryon couplings and their masses, as 



recently claimed in Ref. |[20|| . We find significant reduction in the scaling constant and 
possible deviations from the cubic law. More studies are needed to clarify this important 
issue. 
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